arXiv: 1501.06447vl [astro-ph.IM] 26 Jan 2015 


Mon. Not. R. Astron. Soc. 000, [H fT^ (2015) Printed 27 January 2015 (MN style file v2.2) 


Generalized Formalisms of the Radio Interferometer 
Measurement Equation 

D. C. Price^* and O. M. Smirnov^’^ 

^Harvard-Smithsonian Center for Astrophysics, MS42, 60 Garden Street, Cambridge MA, 01243 United States 
^Department of Physics and Electronics, Rhodes University, PC Box 94, Grahamstown, 6 I 4 O, South Africa 
^SKA South Africa, 3rd Floor, The Park, Park Road, Pinelands, 7405, South Africa 


Accepted 2015 January 12. Received 2014 November 04 


ABSTRACT 

The Radio Interferometer Measurement Equation (RIME) is a matrix-based mathe¬ 
matical model that describes the response of a radio interferometer. The Jones calculus 
it employs is not suitable for describing the analogue components of a telescope. This 
is because it does not consider the effect of impedance mismatches between compo¬ 
nents. This paper aims to highlight the limitations of Jones calculus, and suggests 
some alternative methods that are more applicable. We reformulate the RIME with 
a different basis that includes magnetic and mixed coherency statistics. We present 
a microwave network inspired 2N-port version of the RIME, and a tensor formalism 
based upon the electromagnetic tensor from special relativity. We elucidate the limita¬ 
tions of the Jones-matrix-based RIME for describing analogue components. We show 
how measured scattering parameters of analogue components can be used in a 2N-port 
version of the RIME. In addition, we show how motion at relativistic speed affects the 
observed flux. We present reformulations of the RIME that correctly account for mag¬ 
netic field coherency. These reformulations extend the standard formulation, highlight 
its limitations, and may have applications in space-based interferometry and precise 
absolute calibration experiments. 

Key words: Methods: data analysis — Techniques: interferometric — Techniques: 
polarimetric 


1 INTRODUCTION 


Coherency in electromagnetic fields is a fundamental topic 
within optics. Its importance in fields such as radio as¬ 
tronomy can not be overstated: interferometry and syn¬ 
thesis imaging techniques rely heavily upon coherency the¬ 


ory (Taylor et al. 1999 Mandel & Wolf 1995 Thompson 


|et al. 20041. Of particular importance to radio astronomy 


is the Van-Cittert-Zernicke theorem (vC-Z, Zernicke 19381 


and the radio interferometer Measurement Equation (RIME, 
|Hamaker et al.|[l996 K The vC-Z relates the brightness of a 
source to its mutual coherency as measured by an interfer¬ 
ometer, and the RIME provides a polarimetric framework 
to calibrate out corruptions caused along the signal’s path. 

While the vC-Z theorem dates back to 1938, more re¬ 


cent work such as that of Carozzi & Woan (20091 extends its 


applicability to polarized measurements over wide fields of 
view. The RIME has a much shorter history: it was not for¬ 
mulated until 1996 ( Hamaker et al.[1996 1. Before the RIME, 
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calibration was conducted in an ad-hoc manner, with each 
polarization essentially treated separately. The framework 


was expounded in a series of follow-up papers (Sault et al. 


1996 Hamaker fc Bregman|19^ Hamaker|2006 20061; re¬ 

cent work by Smirnov extends the formalism to the full 
sky case, and reformulates the RIME using tensor algebra 
I Smirno^|201 la|b I. 


This article introduces two reformulations of the RIME, 
that extend its applicability and demonstrates limitations 
with the Jones-matrix-based formalism. Both of these refor¬ 
mulations consider full electromagnetic coherency statistics 
(i.e. electric, magnetic and mixed coherencies). The first is 
inspired by transmission matrix methods from microwave 
engineering. We show that this reformulation better ac¬ 
counts for the changes in impedance between analogue com¬ 
ponents within a radio telescope. The second reformulation 
is a relativistic-aware formulation of the RIME, starting with 
the electromagnetic field tensor. This formalism allows for 
relativistic motion to be treated as instrumental effect and 
incorporated into the RIME. 


This article is organized as follows. In Section 2, we re- 


© 2015 RAS 














































2 D.C. Price and O. Smirnov 


view existing formalisms of the RIME and methods from mi¬ 
crowave engineering. Section 3 defines coherency matrices, 
which are used in Section 4 to formulate general coherency 
relationships for radio astronomy. In Section 5, we intro¬ 
duce a tensor formulation of the RIME based upon the elec¬ 
tromagnetic tensor from special relativity. Discussion and 
example applications are given in Section 6; concluding re¬ 
marks are given in Section 7. 


2 JONES AND MUELLER RIME 
FORMULATIONS 


Before continuing on to derive a more general relationship 
between the two-point coherency matrix and the voltage- 
current coherency, we would like to give a brief overview and 
derivation of the radio interferometer Measurement Equa¬ 
tion of |Hamaker et aT] ( |1996[ ). Our motivation behind this 
is to highlight that Hamaker et. al.’s RIME is a special case 
of a more general (and thus less limited) coherency relation¬ 
ship. 


In their seminal Measurement Equation (ME) paper. 


Hamaker et al. (19961 showed that Mueller and Jones cal¬ 


culuses provide a good framework for modelling radio inter¬ 
ferometers. In optics, Jones and Mueller matrices are used 
to model the transmission of light through optical elements 
I Jones|[T94T Mueller|[l948 1. Mueller matrices are 4x4 ma¬ 
trices that act upon the Stokes vector 


s = {I Q U Vf, (1) 


whereas Jones matrices are only 2 x 2 in dimension and act 
upon the ‘Jones vector’: the electric field vector in a coor¬ 
dinate system such that z-axis is aligned with the Poynting 
vector 

e{r,t) = {e^{r,t) ey{r,t))'^ . (2) 

Jones calculus dictates that along a signal’s path, any (lin¬ 
ear) transformation can be represented with a Jones matrix, 

J: 


eout{r,t) ^ Jein{r,t) (3) 

A useful property of Jones calculus is that multiple effects 
along a signal’s path of propagation correspond to repeated 
multiplications: 

eout{r,t) = Jn-■ ■ J2Jiein{r,t), (4) 

which can be collapsed into a single matrix when required. 

The RIME uses Jones matrices to model the various cor¬ 
ruptions and effects during a signal’s journey from a source 
right though to the correlator. A block diagram for a (simpli¬ 
fied) two-element interferometer is shown in Figure]^ From 
left to right, the figure shows the journey of a signal from a 
source right through to the correlator. The radiation from 
the source is picked up by two antennas, which we have de¬ 
noted with subscript p and q. The radiation follows a unique 
path to both of these antennas; each antenna also has asso¬ 
ciated with it a unique chain of analogue components that 
amplify and hlter the signal to prepare it for correlation. 
Each of these subscripted boxes may be represented by a 
Jones matrix; alternatively an overall Jones matrix can be 
formed for the p and q branches (the dashed areas). 



Figure 1. Block diagram showing a simple model of an inter¬ 
ferometer that can be modelled with the RIME. Radiation from 
a source propagates through free space to two telescopes, p and 
q. After passing through the telescope’s analogue chain, the two 
signals are interfered in a cross-correlator. 


2.1 Hamaker’s RIME derivation 


The derivation of the RIME is remarkably simple and el¬ 
egant. For a single point source of radiation, the voltages 
induced at the terminals of a pair of antennas, p and q are 


Vp{t) = Jpeo{t) (5) 

Vq{t) = Jqeo{t) ( 6 ) 

In its simplest form, the RIME is formed by taking the outer 
product of these two relationships. Note that in their orig¬ 
inal paper, the authors use the Kronecker incarnation of 
the outer product, which we will denote with *. We reserve 
the symbol ® for the matrix outer product of two matrices 
A® B — AB^, where H denotes the Hermitian conjugate 
transpos^ Using the Kronecker outer product, the RIME 
is given by 

(Up ^ p J q) (CQ ^ ^o) — (^ p J q) ^00 (J) 


where Jp and J, are the Jones matrices representing all 
transformations along the two signal paths, and {Jp * Jq) 
is a 4 X 4 matrix. Here, eoo is the sky brightness of a single 
point source of radiation. For a multi-element interferome¬ 
ter, every antenna has its own unique Jones matrix, and a 
RIME may be written for every pair of antennas. 

Due to their choice of outer product, Hamaker et. al. 
arrive at a coherency vector 


epq{rp,rq,T) = {ep(rp,t) * eq{rq,t + t)) 



as opposed to the coherency matrix of | Wolf| ( |1954[ ); this is 
introduced in |3.1| below. The column vector of a point source 
at ro is then eoo; that is, p — q and r = 0. The vector eoo 
is related to the Stokes vector by the transfornj^ 
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The quantity (Jp * Jq) in Eq.[^can therefore be viewed as a 


^ For a discussion on the subtleties of outer product definition, 
see|Smirnov|2011b| A.6.1 

^ Here, j = \/—l, to avoid confusion with current, i, used later. 
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Generalized Formalisms of the RIME 3 


Mueller matrix. That is, Eq. [^can be considered a Mueller- 
calculus-based ME for a radio interferometer. To summarize, 
Hamaker et. al. showed that: 

• Jones matrices can be used to model the propagation 
of a signal from a radiation source through to the voltage at 
the terminal of an antenna. 

• A Mueller matrix can be formed from the Jones terms of 
a pair of antennas, which then relates the measured voltage 
coherency of that pair to a source’s brightness. 

Showing that these calculuses were applicable and indeed 
useful for modelling and calibrating radio interferometers 
was an important step forward in radio polarimetry. 


2.2 The 2x2 RIME 


In a later paper, Hamaker (2000 \ presents a modified formu¬ 
lation of the RIME, where instead of forming the coherency 
vector from the Kronecker outer product (*), the coherency 
matrix is formed from the matrix outer product ((g)): 


Epq — (Cp (g Gq) — 


( 10 ) 


The resulting coherency matrix is then shown to be related 
to the Stokes parameters by Eoo = B , where 


B = 


I + Q U + jV 
U-jV I-Q 


The equivalent to the RIME of Eq. is 

('ypg'Uq) — (S) IJ q — 2Tp(eo: 

or more simply. 


eo) Jg 


Vpq ^JpBjf 


( 11 ) 

H 

q 

( 12 ) 

(13) 


This approach avoids the need to use 4x4 Mueller matri¬ 
ces, so is both simpler and computationally advantageous. 
This form is also cleaner in appearance, as fewer indices are 
required. 

I Smirnov] ( |2011a[ ) takes the 2x2 version of the RIME as 
a starting point and extends the RIME to a full sky case. By 
treating the sky as a brightness distribution B(ct), where a 
is a direction vector, each antenna has a Jones term JpG) 
describing the signal path for a given direction. The visibility 
matrix Vpg is then found by integrating over the entire sky: 

Vpg = f Jp(a)B{a)jlf{a)dn. (14) 

J 471 

This is a more general form of Zernicke’s propagation law. 
Smirnov goes on to derive the Van-Cittert Zernicke theorem 
from this result; we return to vC-Z later in this article. 


2.3 A generalized tensor RIME 


In Smirnov (2011bI, a generalized tensor formalism of the 
RIME is presented. The coherency of two voltages is once 
again defined via the outer product e*ej , giving a (l,l)-type 
tensor expression: 


v^; = FlBpl 


(15) 


This formalism is better capable of describing mutual cou¬ 
pling between antennas, beamforming, and wide field po¬ 
larimetry. In this paper, we focus on a matrix based formal¬ 
ism which considers the propagation of the magnetic field in 
addition to the electric field. We then show that this for¬ 
mulation is equivalent to the tensor formalism presented in 


Smirnov (201 lb I, but is instead in the vector space C°. 


2.4 Microwave engineering transmission matrix 
methods 

All formulations of the RIME to date — including the tensor 
formulation — do not consider the propagation of the mag¬ 
netic field. In free space, magnetic field coherency can be 
easily derived from the electric field coherency. However, at 
the boundary between two media, the magnetic field must be 
considered. Here, we introduce some results from microwave 
engineering which contrast with the Jones formalism. 

In circuit theory, the well-known impedance relation, 
V = ZI relates current and voltage over a terminal pair (or 
‘port’). However, this relation is specific to a 1-port network; 
for microwave networks with more than one port, the matrix 
form [V] = [Z] [7] must be used: 


/Ui\ 

V2 


f Z\l Z\2 

Z 21 

ZlN ^ 


(ii\ 

*2 

\vnJ 


\Zn\ . 

Zmn/ 


\*JV/ 


where Zab is the port-to-port impedance from port a to 
port b, Vn is the voltage on port n, and in is the current. A 
common example of a 2-port network is a coaxial cable, and 
a common 3-port network is the Wilkinson power divider. 

The analogue components of most radio telescopes can 
be considered 2-port networks. The 2-port transmission, or 
ABCD matrix, relates the voltages and currents of a 2-port 
network: 




(17) 


this is shown in Figure If two 2-port networks are con¬ 
nected in cascade (see Figure]^, then the output is equal 
to the product of the transmission matrices representing the 
individual components: 


f ui\ _ /Al 73i\ /A 2 B2\ /1'2\ 

^ VCi DG\C2 772^*2/’ 


(18) 


as is shown in texts such as Pozar (20051. The elements in 
the 2-port transmission matrix are related to port-to-port 
impedances b^j^ 


A = Z\\IZ2\ 

Z 11 Z 22 — ^ 12^21 
^21 

C = I/Z 21 
D = Z22lZ2\. 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 


^ Note that Z 21 = 0 is zero impedance, which is never satisfied 
in real components, and Z 21 = 00 represents an open circuit 
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4 D.C. Price and O. Smirnov 



Figure 2. Top: The ABCD matrix for a 2-port network. In this 
diagram, voltage is denoted with V, and current with I. Bottom: 
Connecting two components in cascade. Diagram adapted from 
|Pozar| l |2005| l 


Like the Jones matrix, the ABCD matrix allows a sig¬ 
nal’s path to be modelled through multiplication of ma¬ 
trices representing discrete components. While the Jones 
matrix acts upon a pair of orthogonal electric field com¬ 
ponents, the ABCD matrix acts upon a voltage-current 
pair at a single port. As Jones matrices do not consider 
changes in impedance (free space impedance is implicitly 
assumed), it is not suitable for describing analogue compo¬ 
nents. Conversely, the 2x2 ABCD matrix cannot model 
cross-polarization response of a telescope. In the section 
that follows, we derive a more general coherency relation¬ 
ship which weds the advantages of both approaches. 


and [rq,t -|- t ), respectively. We may arrange these into a 
single 6x6 matrix Bp, that is equivalent to the time averaged 
outer product of the electromagnetic held column vectors at 
spacetime points (rp,t) and (r,,t-|-r): 


Bp, = 






( 


Mp,\ 

H,J 


(29) 


This matrix fully describes the coherency properties of an 
electromagnetic held at two points in spacetime. We will 
refer to this matrix as the two point coherency matrix. It is 
worth noting that: 


• When Tp = r, and r = 0 we retrieve what [Bergman fc 


Carozzi (20081 refer to as the ‘EM sixtor matrix’. Bergman 


& Carozzi (20081 show this sixtor matrix is related to what 


they refer to as ‘canonical electromagnetic observables’: a 
unique set of Stokes-like parameters that are irreducible un¬ 
der Lorentz transformations. These are used in the the anal¬ 
ysis of electromagnetic held data from spacecraft. 

• For monochromatic plane waves, when Vp = r, and r = 
0, and we choose a coordinate system with 2 in the direction 
of propagation (i.e. along the Poynting vector), Epq becomes 
what Smirnov (2011a I refers to as the brightness matrix, B. 

From here forward, we drop the subscript B = Boo and 
shall refer to this as the brightness coherency to highlight 
its relationship with B. 


3 COHERENCY IN RADIO ASTRONOMY 


We now turn our attention to formulating a more general 
RIME that is valid in a larger range of cases. In this section, 
we introduce the coherency matrices of Wolf (19541, along 
with voltage-current coherency matrices. The following sec¬ 
tion then formulates relationships between source brightness 
and measured coherency based upon these matrices. 


3.1 Electromagnetic coherency 

To begin, we introduce the coherency matrices of [Wolfl 
(19541, that fully describe the coherency statistics of an elec¬ 
tromagnetic field. We may start by introducing e{r,t) and 
h{r,t) as the complex analytic representations of the electric 
and magnetic held vectors at a spacetime point {r,t): 


3.2 Voltage and current coherency 

A radio telescope converts a free space electromagnetic held 
into a time varying voltage, which we then measure after sig¬ 
nal conditioning (e.g. amplihcation and hltering). As such, 
radio interferometers measure coherency statistics between 
time varying voltages. 

One may model the analogue components of a telescope 
as a 6-port network, with three inputs ports and three out¬ 
put ports. We propose this so that there is an input-output 
pair of ports for each of the orthogonal components of the 
electromagnetic held. We can then dehne a set of voltages 
v(t) and currents i{t) 


v{t) = (vx{t) Vy{t) Vz{t)Y (30) 

i(t) = [ixit) iy{t) iz{t)Y. (31) 


e{r,t) = [e^{r,t) ey{r,t) e^{r,t)Y (23) 

h(r,t) = {h,,{r,t) hy{r,t) h,,{r,t))'^ (24) 

The coherency matrices are then dehned by the formulae 
Epq{rp,rq,T) = {{ek(rp,t)ei{rq,t + T))) (25) 

Hpq{rp, rq,r} = {{hk{rp,t)h* *i {rq,t + r))) (26) 

Mpq(rp,rq,T) = {(ek{rp,t)h*{rq,t + T))) (27) 

Npq{rp, r,, r) = {{hk{rp, t)ei (r„ t A r))) . (28) 

Here, k and I are indices representing the x, y, z subscripts 
from Cartesian coordinates. Epq and Hpq are called the elec¬ 
tric and the magnetic coherency matrices, and Mp, and 
Npq are called the mixed coherency matrices. The sub¬ 
scripts p and q correspond to the spacetime points {rp,t) 


In practice, most telescopes are single or dual polarization, 
so only the x and y components are sampled. Nonetheless, it 
is possible to sample all three components with three orthog¬ 
onal antenna elements (Bergman et al. 20051. The voltage 
response of an antenna is linearly related to the electromag¬ 
netic held strength ( jHamaker et al.||1996| ), and the current 
is linearly related to voltage by Ohm’s law, so we may write 
a general linear relationship 


("(^)^ ^ M (e{r,t)\ 

\i{t)) \C d) \h{r,t)) 


(32) 


where A, B, C and D are block matrices forming an over¬ 
all transmission matrix T'. We can now dehne a matrix 
of voltage-current coherency statistics that consists of the 
block matrices 


© 2015 RAS, MNRAS 000, [iHT^ 





















































Generalized Formalisms of the RIME 5 


VpqG) = {{Vk{t)vt{t + T))) 

(33) 

WpqG) = {{ik{t)i*i{t + T))) 

(34) 

KpqG) = + t))) 

(35) 

LpqG) = {{ik{t)vt{t + T))) , 

(36) 


these are analogous to (and related to) the electromagnetic 
coherency matrices abov^ In a similar manner to the two- 
point coherency matrix, we dehne Ypq 



which we will refer to as the voltage-current coherency ma¬ 
trix. This is analogous to the ‘visibility matrix’, Vp,, of 


Smirnov (2011al. 


4 TWO POINT COHERENCY 
RELATIONSHIPS 

Now we have introduced the two-point coherency matrix 
Bpq and the voltage-current coherency matrix Ypq, we can 
formulate relationships between the two. In this section, we 
hrst formulate a general coherency relationship describing 
propagation from a source of electromagnetic radiation to 
two spacetime coordinates. We then show that this relation¬ 
ship underlies both the RIME and the vC-Z relationship. 


4.1 A general two point coherency relationship 

Suppose we have two sensors, located at points rp and r,, 
which fully measure all components of the electromagnetic 
held. Assuming linearity, the propagation of an electromag¬ 
netic held /q = (e(ro,t) /i(ro,t))^ from a point ro to rp 
and rq can be encoded into a 6x6 matrices, Tp and T^: 


fp = Tp/o (38) 

fq = T,/o (39) 

The coherency between the two signals f ^ and f ^ is then 
given by the matrix Bp^: 

~ (/p ®/<j) (40) 

= ((Tp/o) ® (T,/o)) (41) 

= (Tp(/o«)/o)Tf) (42) 

= TpBTf (43) 


we can write this in terms of block matrices 

/ Epq Mpq\ _ ( Ap Bp\ ( Eoo Moo\ (Aq Bq\ 
\Npq Hpq)^\Cp Dp)\Nao H^^)\Cq DqJ 

(44) 

This is the most general form that relates the coherency at 
two points Tp and to the electromagnetic energy density 
at point cq. 

In radio astronomy, antennas are used as sensors to 


^ Spatial location r is no longer relevant as the voltage propa¬ 
gates through analogue components with clearly defined inputs 
and outputs. 


measure the electromagnetic held. Following from Eq. |44[ 
we may write an equation relating voltage and current co¬ 
herency: 

Vp, = Tp (TpBTf )t'". (45) 

As the T' and T matrices are both 6 x 6 , we can are both 
collapse these matrices into one overall matrix. Eq. |45| then 
becomes 


Ypq = TpBTf, (46) 

which is the general form that relates the voltage-coherency 
matrix Ypq to the brightness coherency B. 

Equation]^ is a central result of this paper. It is a gen¬ 
eral case which relates the EM held at a given point in space- 
time to the voltage and current coherencies in between pairs 
of telescopes. In the sections that follow, we show that gen¬ 
eralized versions of the Van-Cittert-Zernicke theorem and 
RIME may be formulated based upon this coherency rela¬ 
tionship, and that the common formulations can be derived 
from these general results. 


4.2 The Radio Interferometer Measurement 
Equation 


By comparing Eq. with Eq. |44| it is apparent that the 
Jones formulation of the RIME is retrieved by setting all 
but the top left block matrices to zero, such that we have 

Ypp = ApEooAf. (47) 


But under what assumptions may we ignore the other entries 
of Eq. |44| ? To answer this, we may note that monochromatic 
plane waves in free space have E and H are in phase and 
mutually perpendicular: 


e{r,t) = {eRr,t) ey{r,t)Y 

h{r,t) = — {-ey{r,t) eRr,t))^ (48) 

Co 


Where cq is the magnitude of the speed of light. In such a 
case, all coherency statistics can be derived from the 2x2 
brightness matrix B. Carozzi & Woan (20091 show that the 
field coherencies can be written 


B = 


( ^pq 

\Npq 


Mpq\ 

Hpq) 


where F is the matrix 



F = 


1 

Co 






(49) 


(50) 


Under these conditions, the rank of B is 2, so the relationship 
in Eq. I^is over constrained. It follows that the 2x2 RIME 
is perfectly acceptable — and indeed preferable to Eq. 

— for describing coherency of plane waves that propagate 
through free space. 

There are numerous situations in which we cannot as¬ 
sume that we have monochromatic plane waves. This in¬ 
cludes near field sources where the wavefront is not well ap¬ 
proximated by a plane wave; propagation through ionized 
gas; and situations where we choose not to treat our field as 
a superposition of quasi-monochromatic components. Most 
importantly, the assumptions that underlie the 2x2 RIME 
do not hold within the analogue components of a telescope, 
where the signal does not enjoy free space impedance. 
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4.3 A 2N-port transmission matrix based RIME 

For a dual polarization telescope, a 4-port description (2-in 
2-out) of our analogue system is more appropriate than the 
general 6-port description. Using the result above and 
Eq. |46| we can write a relationship 


V 


pq 


fAp Bp\ ( B BE^ \ (Aq 
lyCp Dp) \FB FBF'^J 



Here, all block matrices have been reduced in dimensions 
from 3 X 3 to 2 X 2. This version of the RIME retains the 
ability to model analogue components, but uses the approx¬ 
imations of the vC-Z to express B in terms of the regular 
2x2 brightness matrix B. The transmission matrices ma¬ 
trices here ar e similar to t he 2N-port transmission matrices 
as defined by Faria (20021 ^ 

The transmission matrices may be broken down into a 
chain of cascaded components. That is, for a cascade of n 
components, we may write the overall transmission matrix 
as a product of matrices representing the individual compo¬ 
nents: 


( Ap B p\ I A Trip B iip\ f A\p Rip\ 

Cp Dp)^\C^p D^p)"'\Cqp Dip) 

This is similar to, but more general than, Jones calculus. 
In the following section we will explore the difference. 


4.4 Limitations of Jones calculus 

Jones calculus essentially asserts two things. Firstly, it as¬ 
serts that the voltage 2-vector at the output of a dual¬ 
polarization system Vp = (vpi,Vp 2 )'^ is linear with respect to 
the EMF 2-vector e = {ex,ey)'^ at the input of the system: 

Vp = J pC, 

where the Jones matrix Jp describes the voltage transmis¬ 
sion properties of the system. The second assertion is that 
for a system composed of n components, the effective Jones 
matrix is a product of the Jones matrices of the components: 

Vp = Jnp • • ■ tjIpG. (^ 4 ) 

With 2N-port transmission matrices, we instead de¬ 
scribe the input of the system by the 4-vector [e, h]'^, which 
for a monochromatic plane wave is equal to 

(h) ^ iei,,ey,-eylco,ei,lcof = , (55) 

and the output of the system is a 4-vector of 2 voltages and 2 
currents (vp, ip)^, which is linear with respect to the input: 



Note that the output voltage is still linear with respect to 
the input e. Indeed, if one is only interested in the voltage, 
the above becomes 


Vp = (Ap -f BpF)e, (57) 


® We note that our definition here is the inverse of that of Faria 
(input and output are swapped). 


i.e. the system has an effective Jones matrix (i.e. a voltage 
transmission matrix) of 

Jp = Ap + BpF. (58) 

However, the Jones formalism breaks down when the sys¬ 
tem is composed of multiple components. For example, with 
2 components, we may naively attempt to apply Jones cal¬ 
culus, and describe the voltage transmission matrix of the 
system as a product of the components’ voltage transmission 
matrices: 

Jp = (-^2p + B2pF){Aip + BipF), (59) 

i.e. 

Vp = (A2pAip J- B2pFA\p J- A2pB\pF -\- B2pFB\pF)e. 

(60) 

This, however, completely neglects the current transmission 
properties. Applying the 2N-port transmission matrix for¬ 
malism, we can see the difference: 



from which we can derive an expression for the voltage vec¬ 
tor: 


Vp = (A2pAip + B2pCip + A2pBipF -f B2pDipF)e, (62) 

which differs from Eq. in the second and fourth term of 
the sum, since in general 

FAip^Cip, FBipj^Dip. (63) 


To summarize, because Jones calculus operates on volt¬ 
ages alone, and ignores impedance matching, we cannot use 
it to accurately represent the voltage response of an analogue 
system by a product of the voltage responses of its individual 
components. The difference is summarized in Eqs. |60[|63| a 
practical example is given in Sect. |6.l] By contrast, the 2N- 
port transmission matrix formalism does allow us to break 
down the overall system response into a product of the com¬ 
ponent responses, by taking both voltages and currents into 
account. 

Since we have now shown the 2x2 form of the RIME 
to be insufficient, the obvious question arises, why have we 
been getting away with using it? Historically, practical ap¬ 
plications of the RIME have tended to follow the formula¬ 
tion of [Noordar^ ( |1996[ ) , rolling the electronic response of 
the overall system (as well as tropospheric phase, etc.) into 
a single ‘G-Jones’ term that is solved for during calibration. 
Under these circumstances, the 2x2 formalism is perfectly 
adequate - it is only when we attempt to model the in¬ 
dividual components of the analogue receiver chain that its 
limitations are exposed. On the other hand, [^arozzi fc Woan| 
(20091 have highlighted the limitations of Jones calculus in 
the wide-held polarization regime. 


5 TENSOR FORMALISMS OF THE RIME 


Up until now, we have presented our 6x6 RIME using ma¬ 
trix notation. We now briehy discuss how the work presented 
here is closely related to the tensor formalism presented in 


Smirnov (2011b I. 


As is discussed in Smirnov (2011bI, the classical Jones 
formulation of the RIME is in the vector space C^. The 
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formulation proposed by Carozzi fc Woan (20091 is instead 


in C^. In contrast, our Eq. 46 can be considered to work in 


C®; that is, our EMF vector has 6 components: 

e = + X] ■ 


(64) 




The coherency of two voltages is once again defined via the 
outer product e*ej , quite remarkably giving a (l,l)-type ten¬ 
sor expression identical Eq. 9 of Smirnov (2011b I: 

= Jl'BpI. (65) 


In the next section, we show an alternative RIME based 
upon a (2,0)-type tensor commonly encountered in special 
relativity. 


5.1 A relativistic RIME 

Another potential reformulation of the RIME involves the 
electromagnetic tensor of special relativity. The classical 
RIME formulation implicitly assumes that both antennas 
measure the EMF in the same inertial reference frame. If this 
is not the case (consider, e.g., space VLBI), then we must 
in principle account for the fact that the observed EMF is 
altered when moving from one reference frame to another, 
and in particular, that the e and h components intermix. In 
special relativity, this can be elegantly formulated in terms 
of the electromagnetic field tensor, which represents the 6 
independent components of the EMF by a (2-0)-type ten¬ 
sor: 


/ 

0 

-Cx/c 

-Cy/C 

-ezfc 


e^/c 

0 

-hx 

hy 


Cyfc 

hx 

0 

hx 

V 

Cz/c 

— hy 

hx 

0 


( 66 ) 


The advantage of this formulation is that the EMF ten¬ 
sor follows standard coordinate transform laws of special rel¬ 
ativity. That is, for a different inertial reference frame given 
by the Lorentz transformation tensor A“/, the EMF tensor 
transforms as: 


(67) 

The measured 2-point coherency between [Ep] and [Ej] can 
be formally defined as the average of the outer product: 


= 2c^([Ep]“'’[E"]^'’), 


( 68 ) 


where represents the conjugate tensor, i.e. [E^]'*''* = . 

A facto r of 2 is introduced for the same reasons as in lSmirno^ 
1 2011a I, and the reason for will be apparent below. Note 


that the indices in the brackets should be treated as labels, 
while those outside the brackets are proper tensor indices. 

Let us now pick a reference frame for the signal (‘frame 
zero’), and designate the EMF tensor in that frame by [Eq], 
or [Eo(®)] to emphasize that this is a function of the four- 
position X — {ct,x). The [Eo(®)] field follows Maxwell’s 
equations; in the case of a monochromatic plane wave prop¬ 
agating along direction z = (l,z), this has a particularly 
simple solution of 


[Fo{x)] = [Eo(:ro)]e-""^ (69) 
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Let us now consider two antennas located at p and q. 
The 2-point coherency measured in frame zero becomes 




-d-fS ^ 2c^{[Fo{p)r^[F"{q)r^) 

= K,[2F{[Fo{xo)r^[F^{xo)r'}]K«, (70) 

where Kp is the complex exponent of Eq. 1691 and is the 


direct equivalent of the A-Jones term of the RIME (Smirnov 


2011a I. The quantity in the square brackets is the equivalent 


of the brightness matrix B, which we’ll call the brightness 
tensor-. 


B 




= = 2c"([Eo(:ro)]“'’[Eo"(*o)]^"). (71) 




Each element of the brightness tensor gives the coherency 
between two components of the EMF observed in the cho¬ 
sen reference frame (‘frame zero’). Nominally, the brightness 
tensor has 4^ = 256 components, but only 36 are unique and 
non-zero (given the 6 components of the EMF). [Carozzi fc 
Bergman] 120061 show that the brightness tensor may be de¬ 
composed into a set of antisymmetric second rank tensors 
(‘sesquilinear-quadratic tensor concomitants’) that are irre¬ 
ducible under Lorentz transformations. The physical inter¬ 
pretation of the 36 unique quantities within the brightness 
matrix is discussed in [Bergman &: Carozzi] ( |2008| ) , with re¬ 
gards to the aforementioned irreducible tensorial set. 

While the brightness tensor has redundancy not present 
in the tensorial set of Carozzi & Bergman (20061, we shall 


continue to use it as a basis to our relativistic RIME for clar¬ 
ity of analogy to the brightness coherency matrix of Eq. |46[ 
and as it leads to a relativistic RIME for which we can 
define transformation tensors analogous to Jones matrices. 
The redundancy can be described by a number of symmetry 
properties of the brightness tensor: it is (a) Hermitian with 
respect to swapping the first and second pair of indices: 


B‘ 


Y^'yS _ ^'ySap 


(72) 


(b) antisymmetric within each index pair (since the EMF 
tensor itself is antisymmetric, i.e. E“^ = —E^“): 


B' 


y.^^5 _ 


= 


'y.-yd _ 


= -B' 


il3S^ 


(73) 


To see the direct analogy to the brightness matrix, consider 
again the case of the monochromatic plane wave propagating 


along z (Eq. 481. The EMF tensor then takes a particularly 
simple form: 


E“p = 


/ 0 


V 0 


€-X 

0 

0 


-Sy 

0 

0 

Cy 


0 


(74) 


0 / 


and the brightness tensor has only 8^ = 64 non-zero compo¬ 
nents, with the additional ‘0-3’ symmetry property: 


,0/375 

a/BOS 


= B 
= B 


3/375 

0:^35 


0:075 

a/B'yO 


= B 
= B 


0/373 


(75) 


Four unique components can be defined in terms of the 
Stokes parameters: 


^0110 =I + Q 
50120 


^0220 ^ J_Q 
^0210 _-y 


(76) 
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01 02 10 13 


01 -/ II 

02 -I 

10 I -I -I 

13 I -I -I 

20 I 

23 I 

31-1 II 

32 -I 


20 


I 


-I 

-I 

I 


23 31 

-I 

I 

I 

I 

-I 

-I 

-I 

I 


32 


-I 


I 

I 

-I 


Table 1. The non-zero components of the brightness tensor 
unpolarized plane wave. Rows correspond to 

columns to 7(5. 


and conversely, 
_B°iio + 5° 


I = 
U = 


^0120 _|_ ^0210 




V = 


j^OllO _ j^0220 


^0120 _ ^0210 
2i ■ 


(77) 


The other non-zero components of the brightness tensor can 
be derived using the Hermitian, antisymmetry and 0-3 sym¬ 
metry properties. Finally, for an unpolarized plane wave, 
only 32 components of the brightness tensor are non-zero 
and equal to ±7. As these will be useful in further calcula¬ 
tions, they are summarized in Table 

So far this has been nothing more than a recasting of 
the RIME using EMF tensors. Consider, however, the case 
where antennas p and q measure the signal in different iner¬ 
tial reference frames. The EMF tensor observed by antenna 
p becomes 


= (78) 

where [Ap]“ is the Lorentz tensor corresponding to the trans¬ 
form between the signal frame and the antenna frame. Since 
the same Lorentz tensor always appears twice in these equa¬ 
tions (due to F being a (2,0)-type tensor), let us designate 

= [Ap]“,[A,]^, (79) 

for compactness. The measured coherency now becomes 

The equivalent of Jones matrices would be (2,2)-type ‘Jones 
tensors’ so a more general formulation of the above 

would be 




DC^'jS _ 


[jpUABor 


lA 


(81) 


where tensor conjugation is dehned as: 

= (82) 

Note that both the K and A terms of Eq. [^can be consid¬ 
ered as special examples of Jones tensors. The K term can 
be explicitly written as a (2,2)-type tensor via two Kronecker 
deltas: 


(83) 

while the A term is a (2,2)-type tensor by definition (Eq.|79[), 
noting that 

[Ap"l“f,, = [Ap]:f,„ (84) 

since the components of any Lorentz transformation tensor 
A“/ are real. 


Finally, let us note the equivalent of the ‘chain rule’ 
for Jones tensors. If the signal chain is represented by a se¬ 
quence of Jones tensors [Jp,n], ..., [Jp,i] (including K terms, 
A terms, and all instrumental and propagation effects), then 
the overall response is given by 




[Jp,i 


(85) 


Equations |80[[M] and |85| above constitute a relativistic refor¬ 
mulation of the RIME (RRIME). The RRIME allows us to 
incorporate relativistic effects into our measurement equa¬ 
tion. That is, we can treat relativistic motion as an ‘instru¬ 
mental’ effect. Some interesting observational consequences, 
and the relation of the RRIME to work from other fields are 
treated in the discussion below. 


6 DISCUSSION 


The formulations presented here highlight the relationship 
between the seemingly disparate fields of microwave net¬ 
working and special relativity to measurements in radio as¬ 
tronomy. This is a remarkable illustration of how these fields 
are intrinsically related by underlying fundamental physics. 

Indeed, it has long been known that Jones and Mueller 
transformation matrices are related to Lorentz transforma¬ 
tions by the Lorentz group. Wiener was aware as early as 
1928 that the the 2x2 coherency matrix could be written 
in terms of the Pauli spin matrices (Wiener 1928| 19301. 
More recently, Baylis et al. (19931 presents a more general 


geometric algebra formalism that unifies Jones and Mueller 
matrices with Stokes parameters and the Poincare sphere. 


Han et al. (1997b I show that the Jones formalism is a repre¬ 


sentation of the six-parameter Lorentz group; further to this 


Han et al. (1997a I show that the Stokes parameters form a 


Minkowskian four-vector, similar to the energy-momentum 
four-vector in special relativity. 

In contrast, the RRIME presented here is novel as it 
fully describes interferometric measurement between two 
points, instead of simply describing polarization states and 
defining a transformation algebra. Similarly, the 2N-port 
RIME specifically shows how microwave networking meth¬ 
ods can be incorporated into a ME. 

What do we gain from using these more general mea¬ 
surement equations in lieu of the simpler 2x2 RIME? For 
most applications, it will suffice to simply be aware of the 
standard RIME’s limitations, and to work in a piecemeal 
fashion. For example, one can quite happily use Jones ma¬ 
trices to describe free space propagation, and microwave net¬ 
work methods to describe discrete analogue components. An 
overall ‘system Jones’ matrix may be derived to describe in¬ 
strumental effects, but this matrix should never be decom¬ 
posed into a Jones chain. 

Similarly, special relativity describes relativistic effects 
through Lorentz transformations acting upon the EMF ten¬ 
sor. We can treat relativistic motion as an instrumental ef¬ 
fect by using the RRIME, or we can apply special relativistic 
corrections separately, as required. The effect of relativistic 
boosts on the Stokes parameters are considered in|Cocke &| 
Holm|([T972l). 


In the subsections that follow, we present some poten¬ 
tial use cases that highlight the usefulness of the 2N-port 
and relativistic RIME formulations. 
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6.1 Modelling real analogue components 

The 2N-port RIME may have practicality in absolute flux 
calibration of radio telescopes. Generally, interferometer 
data is calibrated by assuming that the sky brightness is 
known, or at least approximately known. If we wish to do 
absolute calibration of a telescope without making assump¬ 
tions about the sky, we must instead make sure that we 
accurately model the analogue components of our telescope. 
A particularly ubiquitous method of device characterization 
within microwave engineering is the use of scattering param¬ 
eters; we briefly introduce these below before incorporating 
them into the 2N-port RIME. 


6.1.1 Scattering parameters 

Voltage, current and impedance are somewhat abstract con¬ 
cepts at microwave frequencies, so engineers often use scat¬ 
tering parameters to quantify a device’s characteristics. 
Scattering parameters relate the incident and reflected volt¬ 
age waves on the ports of a microwave network. The scat¬ 
tering matrix, S, is given by 


Al \ 


/5ii 5 i2 • 

Sln\ 


/vt\ 

4- 

^2 

= 

521 



«2 

V^^n/ 


\5nl ■ • • • 

' Snn/ 




where vf is the amplitude of the voltage wave incident on 
port n, and is the amplitude of the voltage wave reflected 
from port n. 

For a dual polarization system (we will label the po¬ 
larization X and y) with negligible crosstalk, we can model 
the analogue chain for each polarization as a discrete 2-port 
network. Assuming that the analogue chains have the same 
number of components (but not that the components are 
identical), the transmission matrix for each pair of compo¬ 
nents is 


/i. 

0 

B^ 


0 

Ay 

0 

By 

a 

0 

b^ 

0 

V 0 

Cy 

0 

DyJ 


where the elements are from the ABCD matrices of the x 
and y polarizations, and are given by the relations 


A = 

B = 
C = 

b = 


1 + S'l 2 S '21 + S 22 — •S'll(l -I- 522 ) 

2S'i2 

„ 1 + 5ll — S'l 2 S '21 -I- S 22 + S'llS '22 

-2ft^- 

1 —1 + 5ll -|- S'l 2 S '21 -I- S 22 — S'llS '22 
Zo 2Si2 

1 -t- Sii -t- S 12 S 21 — S 22 — S 11 S 22 
2S^^ 


( 88 ) 

(89) 

(90) 

(91) 


Here, Zq is the characteristic impedance of the analogue 
chain, which in most telescopes is set to 50 or 75 ohms. We 
have added tildes to the ABCD parameters, as we are using 
the inverse definition to that in Pozar (20051. 

If the system has significant crosstalk between polar¬ 
izations, the analogue chain is more accurately modelled as 
a 4-port network. In this case the relationships are not so 



Figure 3. Block diagram of a three-state switching experiment. 
In addition to the antenna path, an RF load can in series with 
a noise diode can be selected. The noise diode can be turned on 
and off, giving three possible states: antenna, load, and diode -|- 
load. These three states are used for instrumental calibration. 


simple, and the off-diagonal entries of the block matrices of 
T will no longer be zero. 


6.1.2 Scattering matrix example 

We now present a simple illustration of the differences be¬ 
tween assigning a component a scalar value (as is done in the 
Jones formalism), and by modelling it as a 2-port network. 
Consider a component with a scattering matrix 


/0.1Z0° 0.9Z0°\ 

'v0.9Z0° O.IZOV 


(92) 


for a given quasi-monochromatic frequency. Here, values are 
presented in angle notation to emphasize that they are com¬ 
plex valued. In decibels, the and S 22 parameters have a 
magnitude of -10 dB, and the 512 and ^21 parameters have 
a magnitude of about —0.5 dB. Now suppose we have three 
identical copies of this component, and we connect them to¬ 
gether in cascade. If one only considered the forward gain 
(521), one would arrive at an overall 52itot of 

52itot = 52 i52i52i = 0.729 (93) 


In contrast, using the standard microwave engineering meth¬ 
ods, we can form a transmission matrix (using equations |88f 
|91| above), and then convert this back into an overall matrix, 
5tot. By doing this, one finds 


/0.25Z0° 0.75Z0°\ 

'y0.75Z0° 0.25Z0°J ' 


(94) 


The difference becomes more marked as 5ii and S 22 in¬ 
crease; as 5ii and S 22 approach zero, the two methods con¬ 
verge. 

When designing a component, 5ii and S 22 are generally 
optimised to be as small as possible over the operational 
bandwidtlj^ Nevertheless, their affect on the overall system 
is often non-negligible. 


6.1.3 Absolute calibration experiments 

Now we have shown how scattering parameters can be used 
within the 2N-port RIME, we turn our attention to the chal- 

® A notable exception is RF filters, for which Sn is often close to 
unity out-of-band. In such cases 5ii is strongly dependent upon 
frequency. 
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10 D.C. Price and O. Smirnov 


lenges of absolute calibration of radio telescopes. Measure¬ 
ment of the absolute flux of radio sources is fiendishly hard; 
as such, almost all calibrated flux data presented in radio 
astronomy literature are pinned to the flux scale presented 
in 


Baars et al. (19771, either directly or indirectly (Keller- 
mann|2d09 1. Recent experiments, such as the Experiment to 


Detect the Global EoR Step (EDGES, Bowman & Rogers 


2010), and the Large-Aperture Experiment to detect the 


Dark Ages (LEDA, Greenhill & Bernardi 2012 i, are seek¬ 


ing to make precision measurement of the sky temperature 
(i.e. total power of the sky brightness) as a function of fre¬ 
quency. The motivation for this is to detect predicted faint 
(mK) spectral features imprinted on the sky temperature 
due to coupling between the microwave background and neu¬ 
tral Hydrogen in the early universe. For such instruments — 
and other instruments with wide field-of-views — this flux 
‘bootstrapping’ method is not sufficient. 

For experiments such as EDGES a thorough under¬ 
standing of the analogue systems is vital to control the sys- 
tematics that confound calibration. The calibration strategy 


used for EDGES is detailed in Rogers & Bowman (2012), 


and consists of a novel three-state switching and rigorous 
scattering parameter measurements with a Vector Network 
Analyzer (VNA) to remove the bandpass and characterize 
the antenna. 

Our 2N-port RIME allow for scattering parameters to 
be directly incorporated into a measurement equation that 
relates sky brightness to measured voltages at the digitizer. 
This could be used either to (a) infer the scattering parame¬ 
ters of a device given a known sky brightness, or (b) infer the 
sky brightness from precisely measured scattering parame¬ 
ters. The EDGES approach is the latter, but via an ad-hoc 
method without formal use of a ME. 


6 . 1.4 A three-state switching measurement equation 
The three-state switching as described in [Rogers fc Bow-| 


(20121 involves switching between the antenna and a 


reference load and noise diode, and measuring the resultant 
power spectra {Pa, Pl, and Pd, respectively). This is shown 
as a block diagram in Figure The power as measured in 


each state is then given by: 


Pa = Gks Av{Ta + Trx) 

(95) 

Pl = GkBAv{TL + Tx) 

(96) 

Pd = GkgAviTD +Tl + Tx) 

(97) 


where ks is the Boltzmann constant, Td and Tl are the 
diode and reference load noise temperatures, Trx is the re¬ 
ceiver’s noise temperature, G is the system gain, and Av is 
bandwidth. One can recover the antenna temperature Ta by 


Ta = Td^ - ^+Tl, 


Pd-Pl ' 

where the diode and load temperatures. To and Tl, are 
known. Antenna temperature is then related to the sky tem¬ 
perature by 


Tsky = T^(i-|r|^), 


(99) 


This effect is a prime example why one must not use Jones 
matrices to describe analogue components separately; an 
example showing a standing wave present on three-state 
switch calibrated spectrum from a prototype LEDA antenna 
is shown in Figure]^ 

We may instead write the equations above in terms of 
2N-port transmission matrices, and form MEs for the three 


states: 

Va = 

G(TAB,kyTf -bB,x)G" 

(100) 

Vl = 

G(Bi -bB,x)G" 

(101) 

Vd = 

G(TiBDTf +Bi -bB,x)G". 

(102) 


Here, we have replaced the scalar powers Pa,l,d with corre¬ 
sponding voltage-coherency matrices 'Va,l,d, temperatures 
Th.n.sky.rx are replaced with brightness matrices B^.n.sky.rx, 
and Gks is instead represented by a system gain matrix G. 
We have added a transmission matrix for the antenna T^, 
and a transmission matrix T l for the load in series with the 
noise diode. Note that the relation of Eq. |^is now encoded 
into the ME by the matrix Ta, and that we have dropped 
antenna number subscripts as p=q for autocorrelation mea¬ 
surements. 

It is immediately apparent that the cancellations that 
occur in the ratio of Eq. will not in general occur for 
the equivalent ratio R = (Va — V_l)(V_d — Vl)~^. We can 
however retrieve the result of Eq. by treating the two 
polarizations separately, setting G = GI and T_l = I, with 
TA = TA(i-|rp)i. 

Our equations Eq. |100|l02| allow for both polarizations 
to be treated together, which will be important if cross¬ 
polarization terms are non-negligible. Also, we may expand 
Eq. jlOOj to include ionospheric effects. This ability to com¬ 
bine all effects into a single ME may simplify data analysis 
and improve calibration accuracy for such experiments. 


6.2 RRIME and Lorentz boosts 

To illustrate how the RRIME can be used to describe rel¬ 
ativistic effects, let us first consider the ‘simple’ case of a 
relativistically moving source. Suppose we have an unpolar¬ 
ized point source (Table[^, with antennas p and q located in 
the xy plane perpendicular to the direction of propagation 2 
(so the K component becomes unity), both moving parallel 
to the X axis with velocity v. The Lorentz transformation 
tensor from the signal frame to the antenna frame is then 


where 


[A; 


-Via' 


( ^ 

-Pl 

0 

V 0 


-Pl 0 0 \ 

7 0 0 
0 10 ’ 
0 0 1 / 


p = 


1 


1 


(103) 


(104) 


where the reflection coefficient F = Sn. Failure to account 
for reflections (i.e. impedance mismatch) results in an unsat¬ 
isfactory calibration, due to standing waves within the coax¬ 
ial cable that manifest as a sinusoidal ripple on Ta = Ta{v). 


The Lorentz factor 7 is unity 1 at u = 0, and goes to infinity 
with u —^ c. This case can be analyzed without invoking 
coherencies. Consider the EMF, which in the signal frame is 
a plane wave (Eq. 74). In the antenna frame it becomes 
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(a) . (b) 

Figure 4. (a) Data from a LEDA three-state switched antenna during a testing campaign in December, 2013. These correspond to 

the Pa, Pl, ^'Hd Po of Eqns. |95|97| (b) Residuals after applying Eg. |98| to the data shown in Fig. |4a| and removing an antenna model 
power-law fit. The fast varying sinusoid can be attributed standing waves along the ~300 m of coaxial cable connecting the antenna to 
the back-end electronics; the sharp negative spikes are due to radio interference. 




Noting that 


and 


( 0 

Cx 
ley 
\ hex 


-ex -'ycy -hex 
0 P^iey 76 :, 

-hey 0 By 

-■^Bx -By 0 


(105) 



-h ^ 1 



(106) 


(107) 


0 


and 7 “^ + /3^ = 1, we can see that the EMF in the an¬ 
tenna frame is equivalent to a boost of the original EMF by 
7 (also called Dopplsr boost), coupled to a rotation through 
cj) = — cos“^/3 in the xz plane (i.e. the direction of propa¬ 
gation appears to change, also called rBlativistic aborration). 
Both effects are well-understood in the guise of rBlativistic 
beaming, and explain why, for example, some AGNs exhibit 
asymmetric jets (see e.g. Sparks et al.||1992 l. 

From an RRIME point of view, a more interesting case 
arises when one antenna is moving with respect to the other 
(as is the case in space VLBI). Let’s consider antenna p to 
be at rest with respect to the signal frame (so the [Ap]“f^, 
component becomes the equivalent of unity - the product 
of two Kronecker deltas - and thus may be dropped), and 
antenna to be q moving parall el t o the x axis with velocity 
V. The [JqYf, tensor in Eq. 


81 


is then a product of two 


Lorentz tensors of the form oi Fq. |103| In the absence of 
other effects, the measured coherency becomes 


[Vpg 


iap-yS 




(108) 


or writing it out as an explicit sum for two particular ele- 
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ments of interest (and dropping the pq indices): 


yOllO 


(109) 


7<5 


y0220 


(110) 


75 



Each sum above nominally contains 16 terms, but if we as¬ 
sume an unpolarized point source, then from Table(look¬ 
ing up rows ‘ 01 ’ and ‘ 02 ’) we note that only four components 
of B in each sum are non-zero. Combining this with Eq.|103| 


we get: 




T/0110 T/ 2 q1 2\ T 

V = 7(7 - /3 7 ) = 7, 

(111) 

and 




1/0220 ^ 

(112) 


Doing the same sums for and we arrive at zero. 

This implies that our instrument will measure the Stokes 
parameters as 


-^meas 

•^0110 _j_ •^0220 

2 

_ + 7 

2 

(113) 

^^meas 

■^0110 -^0220 

2 

. 71-7 

2 

(114) 

f^meas 

= 0 


(115) 

Vmeas 

= 0. 


(116) 


Since 7 > 1, we measure boost in the total power I, and 
negative apparent Q (that is, linear polarization perpendic¬ 
ular to the direction of motion). The physical meaning of 
this instrumental Q can be understood in terms of the po¬ 
larization aberration discussed by [Carozzi Woan] ( |2009[ ) : 
because the arriving plane wave is aberrated in the frame of 
antenna q (i.e. no longer appears to propagate along the 2 
axis, but along a slightly different direction), it is measured 
as polarized by dipoles that are parallel to xy. 

Note that this formulation does not incorporate the 
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Doppler shift observed by a moving antenna (we quietly as¬ 
sume that the correlator takes care of correcting for this 
when channelizing), or the problems of clock distribution to 
a relativistically moving observing platform. This will have 
to be addressed in a future work. 

In principle, the effects of Doppler boost and relativis¬ 
tic aberration can be described without invoking the full 
RRIME - traditional Jones calculus still suffices. However 
the RRIME provides a unifying framework that allows these 
effects to be incorporated into a single compact interferomet¬ 
ric measurement equation. This is similar to how the origi¬ 
nal RIME formulation of Hamaker et al. 119961 incorporated 
polarimetric effects that were already described previously 
I Sault et al.||199r I in a compact closed form. 

An even more interesting use case arises when the inci¬ 
dent EME can no longer be described by a plane wave. The 
EMF tensor of Eq. then no longer reduces to Eq. |74[ or 
in other words, the e and h components of the EM are no 
longer mutually redundant. Under Lorentz transformations, 
the e and h components then become intermixed - an an¬ 
tenna measuring e in the rest frame will also measure some 
contribution from h. in a moving frame. It is clear that in 
this case, Jones calculus (which only operates on e) can no 
longer apply, and a full RRIME must be invoked. Practi¬ 
cal applications of this will have to wait for near-held space 
VLBI. 


Conceptually, the latter example is very similar to 
the transmission matrix formalism. In a system where the 
voltage and current components intermix, the voltage-only 
Jones formalism no longer applies, and a formalism incor¬ 
porating both components must be invoked. 


7 CONCLUSIONS 

The radio interferometer Measurement Equation provides a 
powerful framework for describing a signal’s journey from 
an astrophysical source to the receiver of a radio telescope. 
Nevertheless, the Jones calculus it employs is in general not 
sufficient to describe the components within the analogue 
chain of a radio telescope. Within the telescope, methods 
from microwave network theory, such as transmission ma¬ 
trices and scattering parameters, are more appropriate. 

Similarly, the Jones formalism is not able to describe 
relativistic effects. We have presented two reformulations 
that account for mixed and magnetic held coherency in a 
way that is not possible with the Jones formalism. These re¬ 
formulations extend the applicability of the RIME to allow 
it to correctly model analogue components, and to describe 
relativistic effects within the RIME. 
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